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Abstract 

We study the error in quadrature rules on a compact manifold, 

<cD{ Zj }V{f}. 

As in the Koksma Hlawka inequality, V {zj} is a sort of discrepancy of the sampling 
points and V{/} is a generalized variation of the function. In particular, we give 
sharp quantitative estimates for quadrature rules of functions in Sobolev classes. 
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1 Introduction 

In what follows, Ai is a smooth compact d dimensional Riemannian manifold with Rie- 
mannian measure dx, normalized so that the total volume of the manifold is 1, and A is 
the Laplace Beltrami operator. This operator is self-adjoint in L 2 (.M), it has a sequence 
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of eigenvalues {A 2 } and an orthonormal complete system of eigenfunctions {(p\(x)}, 
A(fx(x) = \ 2 (p\(x). The eigenvalues, possibly repeated, can be ordered with increas- 
ing modulus. In particular, the first eigenvalue is and the associated eigenfunction is 1. 
An example is the torus T d = R d /Z d with the Laplace operator — d 2 /dx 2 , eigenval- 
ues [An 2 \k\ 2 } keZ d an d eigenfunctions {exp (2irikx)} k£Zd . Another example is the sphere 
S d = {x G M. d+1 , \x\ — 1} with normalized surface measure and with the angular com- 
ponent of the Laplacian in the space R d+1 , eigenvalues {n(n + d — 1)}^ and eigen- 
functions the restriction to the sphere of homogeneous harmonic polynomials in space. 

A classical problem is to approximate an integral J M f(x)dx with Riemann sums 
N" 1 YljLi f ( z j)> or weighted analogues Ylf=i u jf ( z j)> an d wnat follows will be con- 
cerned with the discrepancy between integrals and sums for functions in Sobolev classes 
W a,p (.M) with 1 < p < +oo and a > d/p. The assumption a > d/p guarantees 
the boundedness and continuity of the function / (x), otherwise / (zj) may be not de- 
fined. As a motivation, assume there exists a decomposition of M. into iV disjoint pieces 
M. = Ui U U2 U ••• U Un and these pieces have measures N^ 1 and diameters at most 
cN~ l l d . Choosing a point Zj in each Uj, one obtains the estimate 



N 

j=1 Jm 



f(x)dx 



N 

<£ / \f(*i)-f(*)\*c< su p (l/(y)-/( 

• =1 JUj \y~x\<cN- 1 / d 



x)\} 



In particular, since functions in W a ' p (M) with a > d/p are Holder continuous of 
degree a — d/p, one obtains 



N 

tf" 1 £/(*;)- / f(^)dx 
,=i Jm 



< cN~( a ~ d l p ^ d 



On the other hand, it will be shown that suitable choices of the sampling points {zj} 
improve the exponent l/p — ct/d to — ct/d and this is best possible. More precisely, the 
main results in this paper are the following: 

(1) For every d/2 < a < d/2 + 1 there exists c > Osuchthatif M = U1UU2U...UUN 
is a decomposition of the manifold in disjoint pieces with measure \Uj\ = Uj, then there 
exists a distribution of points {zj}^ =1 with Zj G Uj such that for every function f(x) in 
the Sobolev space W a ' 2 (M), 



N 

^Ujf{zj) - I f{x)dx 
, =1 Jm 



<c max {diameter (Uj) a } 

l<j<N 1 



fa,2 ■ 



(2) Assume that the points {zj}^ =1 and the positive weights {w,} iV =1 give an exact 
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quadrature for all eigenfunctions with eigenvalues X 2 < r 2 , that is 

N 



fa) 

3=1 



ip\(x)dx 



M 



1 if X = 0, 
if < A < r. 



Then for every 1 < p < +oo and a > d/p there exist c > 0, which may depend on 
M.,p, a, but is independent of r, {zj}J =l and {wj}j =1 , such that 



N 

, =1 J M 



f(x)dx 



< cr~ a 



(3) If 1 < p < 



{ z j}f=i cmd positive weights {ujj}^ = 



oo and a > d/p, then there exist c > and sequences of points 
with 



N 

y^ujjfizj) - / f(x)dx 



< cN~ a/d 



(4) For every 1 < p < +oo and a > d/p there exists c > such that for every distri- 
bution of points {zj}^ =l and numbers {ujj}^ =1 there exists a function f{x) in W a ' p {M) 



with 



N 

E 

3=1 



Ujf {Zj) 



f(x)dx 



M 



> cN- a/d 



An explicit example is the following. The torus T d can be partitioned into N 



n 



congruent cubes with sides 1/n and this partition generates the mesh of points {rr l l/) fl 
T d , which gives an exact quadrature at least for all exponentials exp (2%ikx) with \kj\ < 
n. In this case, (1) and (2) give an upper bound for the error in numerical integration of 
the order of N~ a / d . More generally, if a manifold is decomposed into iV disjoint pieces 
M = Ui U U 2 U ... U U N with diameters < cN~ 1/d , then (1) gives the upper bound N~ a/d . 
Moreover, for every r > there are approximately cr d eigenfunctions with eigenvalues 
A 2 < r 2 and one can choose iV < cr d nodes {zj}^ =1 and positive weights {w/} - =1 which 
give an exact quadrature for these eigenfunctions. Then in this case (2) gives the above 
upper bound N~ a l d . Hence (1) and (2) imply (3), and by (4) this latter is optimal. When 
the manifold is a torus or a sphere these results are essentially known, and indeed there is 
a huge literature on this subject. See il24l for deterministic and stochastic error bounds in 
numerical analysis. In particular, (3) and (4) for p = 2 and for spheres are contained in 
fTI , [fT5l and [fT6ll . For Besov spaces on spheres a result slightly more precise than (3) is in 
ifTTl . while a result slightly weaker than (4) for compact two point homogeneous spaces is 
in [|2ni . See also [10] and, for a survey on related results, lfT4ll and [fT9l . Beside the proofs 
of (1), (2), (3), (4), which are contained in the following section, the paper contains also 
a final section with a number of further results and remarks. Among them it is proved 
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that if a quadrature rule gives an optimal error in the Sobolev space W a ' 2 (Ai), then this 
quadrature rule is optimal also in all spaces W 73 ' 2 (Ai) with d/2 < (3 < a. Moreover, it is 
proved that there is a relation between quadrature rules and geometric discrepancy: 

(5) If dv(x) is a probability measure on Ai, then the norm of the measure dv(x) — dx 
as a linear functional on W a ' 2 (Ai) decreases as a increases. Moreover, if the norm of 
dv(x) - dx on W a ' 2 (Ai) is r~ a , 



f(x)du(x) — / f(x)dx 

M J M 



then for every d/2 < (3 < a there exists a constant c which may depend on a, (3, Ai, but 
is independent of r and dv(x), such that 



f(x)dv(x) 



M 



f(x)dx 



M 



< cr 



(6) Assume that for some r > 1 the discrepancy of the probability measure dv(x) with 
respect to the balls {B (y, 5)} with center y and radius 5 satisfies the estimates 

-- d if 5 < 1 jr, 

if5> 1/r. 



dv(x) 



B(y,S) 



/ dx 


4 


IB{y,5) 





r 



Then for every 1 < p < +oo and a > d/p, there exists a constant c, which may 
depend on a and p, but is independent of dv(x) and r, such that 



f(x)du(x) 



M 



f(x)dx 



M 



cr~ a \\f\\ ma , P ifO<a<l, 
< ( cr- 1 log(l + r) ||/|| w „,, if a 
cr l ll/llw°.p if a > I. 



1, 



Observe that while (1) and (2) hold for specific quadrature rules, (5) is a result for 
arbitrary quadratures. Actually, (5) is only one way, from a to < a. The estimate r~ a 
for an a does not necessarily imply the estimate cr _/3 for f3 > a. Moreover, the sets 
{B(y,S)} in (6) are not precisely geodesic balls, but level sets of suitable kernels on 
the manifold. However, for spheres or compact rank one symmetric spaces these sets are 
geodesic balls. In this case the discrepancy of the measure is the spherical cap discrepancy. 
See 01 or ||23l , and for other relations between quadrature and discrepancy on spheres, 
see also [21. Finally, we would like to point out that our paper is (almost) self contained, it 
does not rely on explicit properties of manifolds or special functions, and it may provide 
a unified vision and simple alternative proofs of some known results. 



2 Main results 



The eigenfunction expansions of functions and operators are a basic tool in what follows. 
The system of eigenfunctions {cp\(x)} is orthonormal complete in L?(Ai) and to every 
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square integrable function one can associate a Fourier transform and series, 



•77(A) = f f{y)vx{y)dy, f(x) = Y,?fW<P 

Jm , 



x(x). 



Since the Laplace operator is elliptic, the eigenfunctions are smooth and it is possible 
to extend the definition of Fourier transforms and series to distributions. In particular, 
the Fourier expansions are always convergent, at least in the topology of distributions. 
One can write the discrepancy between integral and Riemann sum as a single integral 
with respect to a measure d/i(x) = Ylj=i u i&zj{z) ~ dx, with S y (x) the Dirac measure 
concentrated at the point y and dx the Riemannian measure, 

f(x)dfi(x). 

M 



N 

J^Ujfizj) - / f(x)dx 



Then the estimate of the error in the numerical integration reduces to the estimate of 
the norm of a linear functional dp(x) on a space of test functions f(x). Some of the results 
which follow will be stated for generic finite signed measures dfi(x), for measures of the 
form dfx(x) = dv(x) — dx with dv(x) a probability measure, and also for atomic prob- 



ability measures dv{x) = £),- =1 c0j5 Zj (x). The following is an easy and straightforward 
extension to compact manifolds and p norms of some abstract results for reproducing 
kernel Hilbert spaces. See e.g. (D, [[13, [fT3l . 

Theorem 2.1. Let {^(A)} be a numeric sequence indexed by the eigenvalues {A 2 }, with 
{■0(A)} and {^(A) -1 } slowly increasing, that is |^(A)| < a (1 + A 2 )°^ 2 and |^(A) _1 | < 
b(l + A 2 )^^ 2 . Let A(x,y) and B(x,y) be distribution kernels with Fourier transforms 
{ip(\)}and {^(A)- 1 }, 

A ( x i V) = Y1 ^( A VaOVa(2/), B(x, y) = ^2ilj(Xy 1 ipx(x)(px(y). 

A A 

Finally, let f(x) be a continuous function and let dfi(x) be a finite measure on Ai. Then, 
If 1 < P, q < +oo and 1/p + 1/q = 1, 



< 



M 



A{x,y)f{y)dy 



M 



f(x)dp(x) 

M 

V ^| 1 /P 

dx\ 

' M 



B(x, y)dfi(x) 



M 



dy 



1/9 



In particular, when p = q = 2 and B(x, y) = B (y, x) and 



B (x, y) B (y, z)dydfx(z) 



then the above inequality reduces to an equality. 
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Proof. The assumptions {^(A)} and {^(A) -1 } slowly increasing simply imply that the 
kernels A{x,y) and B(x,y) are tempered distributions. In what follows the pairing be- 
tween a test function and a distribution is denoted with an integral, even when the distri- 
bution is not a function and the integral is divergent. Let 

A(x,y)f(y)dy = J2^)J r fW^), 

A 

B(x,y)dn(y) = 5>(Ar\FMA)v?A(*). 

M x 

These operators are one the inverse of the other, 



M 



M 

A(x,y) 



B(x,y)[ / A(y,z)f(z)dz)dy 



M 



M 



B(y,z)f(z)dz dy. 



M 



In particular, by Holder inequality with 1/p + 1/q — 1, 



f(x)dfx(x) 



< 



M 



M 



M JM J M 



B(x, y)A(y, z)f(z)dfj,(x)dydz 



J 



A(y,z)f(z)dz 



dy 



i/p 



M 



B(x,y)dfj,(x) 



M 



dy 



Finally, when p = q = 2 the Cauchy inequality reduces to an equality if the functions are 
proportional. Indeed, if B (x, y) = B (y, x) and 



M JM 



f(x) 

then one easily verifies that 

f(x)dfi(x) = 



B (x,y)B (y,z)dydn(z), 



M 



M 



A{x,y)f{y)dy 



M 



M 

1/2 



B(x,y) dfi(x) 



M 



dy, 



dx 



1 M 



B (x,y)dfj,(x) 



<M 



1/2 



dy 



Hence, when p = q = 2 for this function the inequality in the theorem reduces to an 
equality. □ 

In what follows the operators with kernels A(x, y) and B(x, y) will be powers of the 
Laplace Beltrami operator (/ + A) ±a ^ 2 . 

Definition 2.2. The Sobolev space W a ' p (M), — oo < a < +oo and 1 < p < +oo, 
consists of all distributions on M with (/ + A) a ^ 2 f(x) in L p (M), that is with 

i/p 



M 



£(l + A 2 ) a/ V/(A)v^) 



dx 



< +oo. 
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An equivalent definition is the following. 
Definition 2.3. Let B a (x, y), — oo < a < +00, be the Bessel kernel 

B"(x,y) = (1 + \ 2 Y al \x{x)^M. 
x 

A distribution f(x) is in the Sobolev space W a,p (/A) if and only if it is a Bessel potential 
of a function g(x) in L p (/A), 



B a (x,y)g(y)dy. 



M 



Moreover, \\f\\ Wa , P = \\g\\ LP . 
In particular, when p — 2, 



£(i+A 2 n^(A)r 



1/2 



Another equivalent definition is a localization result: A distribution f(x) is in W a ' p (/A) 
if and only if for every smooth function g (x) with support in a local card x = ip{y) : M. d -w 
A'l, the distribution g(ip(y))f(ip(y)) is in W a ' p (lR d ) . In particular, if a is a positive even 
integer, then f(x) is in W a ' p {M) if and only if the p th power of f(x) and of A a ^ 2 f(x) 
are integrable. Moreover, distributions in W a,p (/A) with a > d/p are Holder continuous 
of degree a — d/p. When applied to functions in Sobolev classes, Theorem [27T1 gives the 
following corollary. 

Corollary 2.4. (1) IfB a (x,y) = £ A (1 + X 2 y a/2 <p x (x)<p x (y) is the Bessel kernel, if 
dp(x) is a finite measure on M., and if f(x) is a continuous function in W a ' p (M.), with 
1 < Pi Q < +00 and 1/p + 1/q = 1, then 



f(x)dp(x) 



M 



< 



M 



B a (x, y)dp,(x) 



M 



dy 



1/9 



If a > d/p then the above integrals are well-defined and finite. On the contrary, the spaces 
W a ' p (M) with a < d/p contain unbounded functions and, if the measure dp(x) does not 
vanish on the set where f(x) = 00, then f M f(x)d/i(x) may diverge. 
(2) When p = q = 2 then the above inequality simplifies, 

1/2 

f(x)dp(x) 



M 



< 



B 01 (x, y) dfi(x)d/j(y) 



M JM 



Equivalently, by the Fourier expansion of the Bessel kernel, 

f(x)dfi(x) 



M 



^{E( 1 + a2 )^i^( a )I 2 } 



1/2 
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Moreover, with f(x) = J M B 2a (x, y) dfj,(y) the above inequalities reduce to equalities. 

(3) Ifdfi(x) = dv(x) — dx is the difference between a probability measure dv(x) and 
the Riemannian measure dx, then 

1/2 



f(x)dis(x) 

'M 

Equivalently, 



f(x)dx 



M 



< 



B 2a (x,y) du(x)dv(y) - 1 



M JM 



f(x)dv(x) 



M 



I f(x)dx 


-4 


> M 





<{E( 1 + a2 )^i^( a )I 2 } 

I A>0 J 



1/2 



Proof. (1) is an immediate corollary of Theorem 12. 1[ In order to prove (2), observe that 

/ B a (x,y)BP(y,z)dy = B a +> 3 (x,z). 

JM 

Moreover, this Bessel kernel is real and symmetric. Hence, 



M 



B a (x,y)dfi(x) 



M 



dy 



M JM JM 



B a (x, y)B a (z, y)dyd/i(x)dfj l (z) 



B a (x, z)d/j,(x)d/j,(z). 

>M JM 

(3) is a corollary of (1) and (2). Indeed, since B 2a (x, y) = B 2a (y, x) and j M B 2a (x, y) dy 
1, it follows that 

B 2a (x, y) (dv(x) — dx) (du(y) — dy) 



B 2a (x,y) dv(x)dy 



M JM 

2a 



B 2a (x,y) dv{x)du{y) 



M JM 



M JM 



-[ [ B 2a {x,y)dxdu(y)+ [ ! B 2a (x,y)dxdy 

J M JM JM JM 



B 2a (x,y) dv{x)dv{y) - 1. 



M JM 



Finally, by Sobolev imbedding theorem, functions in W a ' p (Ai) with a > d/p are contin- 
uous and all the above integrals are well-defined and finite. This also follows from Lemma 
and Remark [331 below. □ 

The above corollary leads to estimate the energy integrals 



M 



B a (x, y)dp(x) 



M 



dy 



1/9 



B 2a (x,y) dv{x)dv(y) - 1 



' M JM 



} 1/2 =\j2(i+\r a \^wt 

' l A>0 



1/2 
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By the last formula, the energy attains a minimum if and only if Tv (A) = for 
all A > 0, and this gives the Riemannian measure dx. The meaning of the corollary is 
that measures with low energy are close to the Riemannian measure and they give good 
quadrature rules. In order to give quantitative estimates for the above integrals, one has to 
collect some properties of the Bessel kernels. The norm of the function y ~» B a (x, y) in 
W 7 ' 2 (M) is 

\\B a (x,.)\\ w , 2 = l^2^ + ^T~lMx)\ 2 

By Weyl's estimates on the spectrum of an elliptic operator, see Theorem 17.5.3 in 
lfi~8l . for every r > 1 there are approximately cr d eigenfunctions (p\(x) with eigenvalues 
A 2 < r 2 and J2\< r I^A(a:)| 2 < cr d . It then follows that the norm in W 7 - 2 (M) of B a (x, y) 
is finite provided that 7 < a — d/2 and, by Sobolev imbedding theorem, it also follows 
that B a (x, y) is Holder continuous of degree 5 < a — d. Indeed, we shall see that a bit 
more is true: B a (x, y) is Holder continuous of degree a — d. 

Lemma 2.5. The heat kernel W(t,x,y) = J2\ ex P ( — ^ 2 t) i P\i. x ) L P\{y)' which is the 
fundamental solution to the heat equation d/dt = —A on M + x M., is symmetric real 
and positive, W(t,x,y) = W(t,y,x) > for every x,y G M. and t > 0. Moreover, 
for every m and n there exists c such that, if \x — y\ denotes the Riemannian distance 
between x and y, and V the gradient, 

f \V m W(t,x,y)\ < cH d+m )/ 2 (1 + \x-y\/Vt)~ n if0<t<l, 
\ \V m W{t,x,y)\ < c ifl<t<+oo. 

Proof. All of this is well known. The idea is that heat has essentially a finite speed of 
propagation and diffusion in manifolds is comparable to diffusion in Euclidean spaces, at 
least for small times. The heat kernel in the Euclidean space W d is a Gaussian, 

W(t,x,y)= [ exp(-47r 2 t|e| 2 )exp(27rz(x- 2 /)0^ 

= (47rt)" d/2 exp (- \x - y\ 2 /At) . 

By the Poisson summation formula, the heat kernel on the torus T d = M d /Z d is the 
periodized of the kernel in the space, 

exp (— 47r 2 |A;| 2 t) exp {2nik {x — y)) 

kez d 

= {^ty d/2 exp (- \x - y - k\ 2 /At) . 

kez d 

When x is close to y and t is small, the main contribution to the sum comes from the term 
with k = 0, 

W (t, x, y) w (Anty d/2 exp (- \x - y\ 2 /At) . 
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The remainder gives a bounded contribution, 



(4nty d/2 exp(-\x -y - k\ 2 /At) 

fcez d -{o} 



< c. 



Analogous estimates hold for the derivatives. This proves the lemma for the torus. The 
heat kernel on a compact manifold is similar, in particular it has an asymptotic expansion 
with euclidean main term. See e.g. O, Chapter VI. More precisely, by the Minakshisun- 
daram Pleijel recursion formulas, there exist smooth functions {uk (x, y)} such that, if t 
is small and \x — y \ denotes the distance between x and y, 

n 

W (t, x, y) « (Ant)^ 2 exp (- \x - y \ 2 /At) ^ t k u k (x, y) + O (t n+1 ) . 

k=0 

On the contrary, W (t, x, y) — ► 1 when t — > +00. The estimates on the size of this 
kernel and its derivatives are a consequence of this asymptotic expansion. The positivity 
W (t, x,y) > is a consequence of the maximum principle for heat equation and the 
symmetry W (t, x,y) = W (t, y, x) follows from this positivity and the eigenfunction 
expansion. □ 

Lemma 2.6. (1 ) The Bessel kernel B a (x, y) with a > Ois a superposition of heat kernels 

W(t,x,y): 

r+00 

B a {x, y)=T (a/2)" 1 / t a/2 ~ l exp (-t) W (t, x, y) dt. 

Jo 

(2) The Bessel kernel B a (x, y ) with a > is real and positive for every x, y G M., 
and it is smooth in {x 7^ y}. Moreover, for suitable constants < a < b, 

a \x — y\ a d < B a (x, y) < b \x — y\ a ~ d ifO < a < d, 
a log (l + I a; — y] 1 ) < B a (x, y) < b log (l + \x — y\ X ) if a = d, 
a < B a (x, y) <b if a > d. 

(3) If d < a < d+ 1, then B a (x, y) is Holder continuous of degree a — d, that is there 
exists c such that for every x,y,z G Ai, 

\B a (x,y) - B a (x,z)\ <c\y-z\ a - d . 

(4) If d < a < d + 2, then there exists c such that for every x, y G Ai, 



\B a (x,x) -B a (x,y)\ <c\x-y 



a—d 
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Proof. When the manifold is a torus and the eigenfunctions are exponentials the proof is 
elementary. The Bessel kernel in the torus T d is an even function and it is sum of cosines, 

B a (x, y) = ^2(l + 4vr 2 \k\ 2 ) ~ a/2 exp (2nikx) exp (-2mky) 
kez d 

= {l + 4:Tr 2 \k\ 2 y a/2 cos(2irk(x-y)). 

kez d 

Hence, 

B a (x, x) - B a (x, y) = 2j2( l + ^ 2 |fc| 2 )~ a/2 sin 2 (nk {x - y)) 

kez d 

<2n 2 \x-y\ 2 J2 \k? {I + ^ 2 \k\ 2 Y a ' 2 + 2 ^ (l + 4vr 2 \k\ 2 y a ' 2 

fcl^lx-j/p 1 fc|>|x-j/| _1 
ice— d 

1 2 log 

|2 







\x 


-y\ 






\x 


-y\ 




i: 


\x 


-y\ 



if d < a < d + 2, 



\x 



y\ x ) if a = d + 2, 



Also observe that the series which defines B a (x, x) — B a (x, y) has positive terms and the 
above inequalities can be reversed. This proves (4) for a torus, and the proof of (3) and (2) 
is similar. A proof for a generic manifold follows from the representation of Bessel kernels 
as superposition of heat kernels and the estimates in the previous lemma. In particular, (1) 
follows from the identity for the Gamma function 

/■+oo 

(l + A 2 )~ Q/2 = r(a/2)- 1 / t a / 2 - 1 exp(-t(l + A 2 ))rft. 

Jo 

By Lemma 1131 for every n, 

f cfr-W \ x - y\~ n if < t < \x -y\ 2 , 
< W (t, x, y) < I ct- d ' 2 if \x-y\ 2 < t < 1, 
{ c if t > 1. 

Hence, if < a < d and n > d — a, 

r+oo 

B a {x, y) = T (a/2)" 1 / t a ' 2 ~ l exp (-t) W {t, x, y) dt 

Jo 

f\x — y\ 2 rl f+oo 

<c\x-y\- n & +n -^ 2 - l dt + c / tW^dt + / t a ' 2 - 1 exp (-t) dt 

Jo J\x-y\ 2 Jl 

^ l \a—d 

< c \x — y\ 

Indeed one can easily see that these inequalities can be reversed. Hence B a (x,y) ~ 
c\x — y\ a ~ d . This proves (2) when < a < d, and the proofs of the cases a = d 
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and a > d are similar. Also the proof of (3) is similar. Write 

B a (x,y)-B a (x,z) 

f a/2-l exp (_ f ) fyy ^ x ^y)_ W ^ 2 )) ^. 

ct~ d/2 if < t < |y- z\ 2 , 



= T(a/2)- 1 
Then recall that, by Lemma [231 



\W(t,x,y)-W(t,x,z)\<{ ct-V+W \y - z\ if \y - z\ 2 < t < 1, 

c | y — z| if i > 1. 



Hence, 







+C |j/ — 2! 



5 a (x,y) - B a (x,z)\ < c 
t 

< c\y — z 



t (a-d)/2-l exp df 



1 r+oo 

(a-d-l)/2-l_ , .I, , _ _| / W2-1 



exp (-i) tft + c|y - z\ J t a/2_1 exp (-t) 

a—d 



Finally, the estimate for \B a (x, x) — B a (x, y) | in (4) is analogous to the previous one, but 
it holds in a larger range of a. It suffices to observe that W (t,x, y) is stationary at x = y 
and it satisfies the estimates 

( ct- d / 2 ifO < t < \x-y\ 2 , 
\W (t,x,x) -W (t,x,y)\ < < ct- d / 2 ^ \x - y\ 2 if \x - y\ 2 < t < 1, 

^ c \x — y\ 2 if t > 1. 

□ 

The following is Result (1) in the Introduction. 

Theorem 2.7. For every d/2 < a < d/2 + 1 there exists c > with the following 
property: If M. = U\ U U 2 U ... U Un is a decomposition of M. in disjoint pieces with 
measure \Uf\ = Uj, then there exists a distribution of points {zj}^ =l with Zj E Uj such 
that 



N 

J2^f(zj) - \ 

j= l JM 



f(x)dx 



< c max {diameter (Ui) a \ 

l<j<N 



2 (M) 



Proof By Corollary El(3), with dv(x) = Ef=i (a) 
f(x)dv(x) 



/ f(x)dx 




' M 





■>.)-- 

N N 



1/2 



i=l j=l 
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It suffices to compute the average value of Yli=i ^2j=i ^i^jB 2a fa, Zj) — 1 on U\ x U 2 x 
... x Un with respect to the probability measures uj~ l dzj uniformly distributed on Uj. 
First observe that 



N 



Then, 



Ff^fe 1 / ••• / dz 1 ...dz N = l, 
Kk=1 J JUi Ju N 

1= [ [ B 2a (x, y) dxdy = VV / / B 2a (x, y) dxdy. 

J M JM i=1 j=1 JU t JUj 



y Uj / B 2a (zj, Zj) dzj + 2_. / / B 2a (zi, Zj) dzidzj 
j ■>> , ' ' i+ j Jujvj 

J2 [ I B 2a {x,y)dxdy- VV / / B 2a (x,y)dxdy 
= E/ / {B 2a (x,x)-B 2a (x,y))dxdy. 

JUi JUi 



j = l " u 3 " u 3 

Since, by Lemma[U(4), \B 2a (x,x) - B 2a (x,y)\ < c \x - y\ 2a ~ d when d < 2a < d+2, 
and since = \Uj\ < c diameter (Uj) d , 

N 



V / / \B 2a (x,x)-B 2a (x,y)\dxdy 

j=l JU 3 JUj 

N 

< I Uj | 2 sup { | B 2a (x, x) - B 2a (x, y) | , x, y G Uj } 

3=1 

N N 

— \ Uj\ 2 diameter (Uj) 2a ~ d < \ Uj \ diameter (U 



,2a 

3=1 3=1 

□ 

For the next result we shall need estimates for partial sums of Fourier expansions of 
the Bessel kernels. 

Lemma 2.8. Let x (A) be an even smooth function on — oo < A < +oo with support in 
1/2 < |A| < 2 and let 

P a (r,x,y) = £x(A/r) (l + A 2 p /2 tp x {x]fifij). 
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Then for every n > there exists c such that for every r > 1 and x,y G M, 

\P a (r,x,y)\ < cr d ~ a (I + r \x - y\y n . 

Proof. The numerology behind this estimate is quite simple. The approximation of the 
Bessel kernel B a (x, y) by linear combinations of eigenfunctions with eigenvalues A 2 < r 2 
is localized and only points x and y with \x — y\ < l/r really matter. In particular, since 
B a (x, y) is smooth away from the diagonal, at distance \x — y\ < l/r the approximation 
is rough, but at distance \x — y\ > l/r it is quite good. The analogue of P a (r, x, y) in the 
Euclidean setting is the kernel 

Q(r,x-y)= [ x(2n\^\/r)(l+An 2 \^\ 2 y a/2 ex V (2m(x-y)Od^ 

JR d 

= r d [ X (2n\Z\)(l + 4n 2 r 2 \Z\ 2 y a/2 exp(2mr(x-y)0dC 

JR d 

Since X (2-7T |£|) has support in 1/2 < 2ir |£| < 2, for every r and x,y G M d one has 
* [ x(^\Z\){l + 4n 2 r 2 \Z\ 2 y a/2 exp(2mr(x-y)0dt 

JR d 

<r d - a [ {2-n\i\)- a \x{2-K\i\)\di<cr d - a . 

JR d 

This estimate can be improved in the range \x — y\ > l/r. Indeed, an integration by parts 
gives 

r d [ X (27r|e|)(l + 47r 2 r 2 |e| 2 )^ /2 exp(2«r( a ;-y)0^ 

jR d 

= r d [ X (2tt (1 + 4ttV |e| 2 ) " Q/2 A£ ((4ttV \x - y\ 2 ) _B exp (2mr (x - y) fl) df 
= r d (4vr 2 r 2 |x - y| 2 p / exp (2vur (x - y) A^ ( x (2vr (l + 4vr 2 r 2 |£| 2 p /2 ) df. 



Hence, 



* f x(^\Z\)(l + 4n 2 r 2 \Z\ 2 y a/2 exp(2mr(x-y)0dt 

JR d 

< r d (4n 2 r 2 \x - y\ 2 y n [ ( X (2vr \£\) (l + 4ttV |£| 2 p /2 

JR d V 



Now it suffices to transfer these estimates from the Euclidean space to the manifold. For 
the torus, this can be done via the Poisson summation formula. If Q (r, x — y) is the 
truncated Bessel kernel in R d defined above, then the truncated Bessel kernel in T d is 

X (2vr \k\ /r) (1 + 4tt 2 \h\ 2 y a/2 exp (2nik (x - y)) = ^ Q (r, x - y + k) . 
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When \xj — yj\ < 1/2, the main term in the last sum is the one with k = 0, while the 
contribution of terms with k ^ is negligible, 

\Q(r,x-y)\ < cr d ~ a (1 + r \x - y\)~ n , 
\Q{r,x-y-k)\<cr d - a - n . 

fcez d -{o} 

Finally, the estimate for the truncated Bessel kernel on a generic manifold can be obtained 
by transference from M. d via pseudodifferential techniques. For more details, see e.g. [|27l 
Chapter XII, or 0. □ 

The following is a result on the homogeneity of measures which appear in quadrature 
rules and it gives sharp estimates of the discrepancy of such measures. Similar estimates 
on spheres are in [0. 

Lemma 2.9. Assume that dv[x) is a probability measure on Ai with the property that 
for every eigenfunction tp\(x) with eigenvalues X 2 < r 2 , 



(p\(x)dv(x) = / (f\(x)dx. 
M Jm 

Then for every n there exists c, which may depend on n and Ai, but is independent of r 
and dv[x), such that for every measurable set Q in A4, 

/ dv(x) — / dx < c / (1 + rdistance {x, dil})~ n dx. 
Jn Jn Jm 

In particular, the discrepancy between the measures du(x) and dx with respect to balls 
{\x — y\ < s} is dominated by 



dv{x) 

{\x-y\<s} 



I dx 




cr 


'{\ x ~y\< s } 




cr 



~ d ifs < l/r, 
^s^" 1 ifs > l/r. 



Proof. It is proved in ifTTl that given n, there exists c such that for every measurable set 
in Ai and every r > there exist two linear combinations of eigenfunctions A(x) = 
Xa<r a W V 9 a(^) and B(x) = J2\< r b W V 9 a(^) which approximate the characteristic 
function xn( x ) f rom above and below, 

A(x) < Xn(x) < B(x), B(x) — A(x) < c (1 + rdistance {x, d£l})~ n . 

In particular, the properties of the function A(x) and of the measure du(x) give 



> 



M 



dis(x) > / A(x)dv(x) = / A(x)dx 
Jm Jm 

Xn(x)dx — c / (1 + rdistance {x, d£l})~ n dx. 
Jm 
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Similarly, by the properties of B(x) and du(x), 



du(x) < / B(x)dv(x) = / B(x)dx 
<> Jm Jm 



< / xn(x)dx + c / (1 + rdistance {x, d£l}) n dx. 
'm Jm 



□ 



Lemma 2.10. Assume that du(x) is a probability measure on Ai which gives an exact 
quadrature for all eigenfunctions <f\(x) with eigenvalues A 2 < r 2 , 

/ (p\(x)du(x) = / (p\(x)dx. 
Jm Jm 

Ifl < q < +oo and a > d (1 — 1/q), then there exists c, which may depend on q, a, M., 
but is independent of r and dv{x), such that 



M 



B a {x,y)dv{x) - 1 



M 



dy }> < cr~ a . 



Proof Let x (A) be an even smooth function on — oo < A < +oo with support in 1/2 < 
|A| < 2 with the property that ZlJ^oo X (2 _j A) = 1 for every A ^ 0. Also, let 

P a (s, x,y) = J2x (A/s) (1 + A 2 ) - a/2 



Hence, B a (x,y) = 1 + Ylj^-oo P a {^ i x i y)- Since du(x) annihilates all eigenfunctions 
with < A < r, it also annihilates all P a (2 J , x, y) with 2 J < r/2 and this gives 



B a {x,y)du(x) - 1 



M 



f I J- P a (2^x,y) \ du{x). 

JM \v>r/2 ) 



When q = 1, by Lemma [2781 



M 



P Q (s, x, y)dv(x) 



M 



dy 



<cs d ~ a / {l + s\x-y\)~ n dv{x)dy 
Jm Jm 

< cs~ a sup if s d (1 + s \x - y\)~ n dy] < cs~ a . 

x&M I J M J 
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When q = +00 and s > r and n > d, by Lemma Upl and Lemma |2~9l 



sup 

yeM 

< cs d ~ a sup 



P a (s, x, y)dv{x) 



M 



yeM Um 

< cs d ~ a sup 



(1 + s \x — y\) n dv(x) 



dv(x) 



y£M \J {\x-y\<l/r} 



{+00 „ 
V (2?s/r)~ n / 
i=0 ' H\*-v 



dv(x) 



-y\<2i/r} 

< cs da r^ d + cs d ~ a ~ n r n ~ d < cs d ~ a r~ d 



Hence, when s > r and 1 < q < +00, by interpolation between 1 and +00, 



M 



P a (s, x, y)du(x) 



M 



dy 



< sup 

yeM 



P a (s, x, y)du(x) 



M 



(9-l)/9 



M 



P a (s, x, y)dv(x) 



M 



dy 



1/9 



< cs d(l-l/g)-o r -d(l-l/g)_ 



When a > d(l — 1/q) these estimates sum to 

/ / B a (x,y)du(x) -1 



< 



A1 



P a (2 j ,x,y)du(x) 



M 



E 

2J >r/2 

< C r" d(1-1/9) 2 i(d(1_1/,2) " a) < 

2J>r/2 



1/9 



□ 



Finally, the existence of exact quadrature rules associated to any system of continuous 
functions is a simple result in functional analysis, or in convex geometry. See Theorem 
3.1.1 in [26], or [25], or [8] for explicit constructions on spheres. 

Lemma 2.11. Given any number (pi(x), (p2(x),..., <f n ( x ) of continuous functions on A4, 
there exist points {Zj}. =1 in M. and positive weights {wj}- =1 with J2j =1 u)j = 1> sucn 
that for every ipi(x), 

N 



/ (p^dx = fa) ■ 

Jm , =1 



If the functions <fj(x) are real one can choose N < n + 1, and if these functions are 
complex one can choose N < 2n + 1. 
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Proof. Define 



E = I ${x)dx 

i M 



= (ipi(x),ip 2 (x),...,ip n (x)) 

(fi(x)dx, / cp 2 {x)dx, / Lp n (x)dx 



, . . . , 

M JM JM 



If all functions (fi(x) are real valued, then and E are vectors in W 1 . If the (fi(x) are 
complex, then and i£ can be seen as vectors in IR 2r \ Moreover, E is in the convex 
hull of the vectors with x G .M. It then follows from Caratheodory's theorem that 
E is also a convex combination of at most n + 1 vectors in the real case, or 2n + 1 
in the complex case, £ = J2f=i u j® i z j) w i m w i > an d SjLi w j = 1- ^ 

The following is Result (2) in the Introduction. 

Theorem 2.12. Assume that the probability measure dv[x) on M. gives an exact quadra- 
ture for all eigenfunctions (p\{x) with eigenvalues A 2 < r 2 , 

(fx(x)du(x) = / (f\(x)dx. 

M JM 

Then, for some constant c independent of dv(x) and r and for every function f(x) in 
W Q ' P (M) with 1 < p < +oo and a > d/p, 



f(x)di/(x) — / f(x)dx 
M Jm 



< cr 



Proof. By Corollary 12.41 (1) with djj,(x) = dv(x) — dx, 



f{x)djj,{x) 



M 



< 



M 



B a (x,y)du(x) 



M 



dy 



1/9 



By the assumption J M Lp\(x)dfi(x) = for every A < r, and Lemma |2.1Q[ 



M 



B a (x,y)du(x) - 1 



M 



dy 



< cr 



□ 



The above theorem has as corollary Result (3) in the Introduction. 



Corollary 2.13. Ifl<p< +oo and a > d/p, then there exists c > with the property 
that for every N there exist sequences of points {zj } =1 and non negative weights {uj }f =v 
such that for every function f(x) in W a ' p (M), 



N 

y^Ujfizj) - / f{x)dx 



< cN~ a/d 
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Proof. By Weyl's estimates on the spectrum of an elliptic operator, see Theorem 17.5.3 
in [fT8l . for a given r there are approximately cr d eigenfunctions ip\(x) with A < r. The 
corollary then follows from Lemma [2.1 H and Theorem 12.121 with r = N l l d . □ 



This corollary for the sphere is contained in nVni . Finally, easy examples show that the 
above estimates for the error in approximate quadrature are, up to constants, best possible. 
Again, see lfT5l for the case of the sphere. In particular, the following is Result (4) in the 
Introduction. 

Theorem 2.14. For every 1 < p < +oo and a > there exists c > with the fol- 
lowing property: For every distribution of points {zj}^ =1 there exists a function f(x) in 
W a,p {M) which vanishes in a neighborhood of these points and satisfies 

, a , p < cN a/d , [ f{x)dx = 1. 



JM 

As a consequence, for every distribution of points {zj}^ =1 and complex weights {ujj}^ =1 , 
there exists a function f(x) with 

N 



yZwjffa) - / f(x)da 

,=1 JM 



> cN~ a/d 



Proof. If e is small, then one can choose 2N disjoint balls in M. with diameters sN^ 1 ^ 
and, given iV points {zj}, at least iV balls have no points inside. On each empty ball 
construct a bump function ipj(x) supported on it with 



||^ll w ^< cAra/d " 1/p > / il>j(x)dx = N- 
Define f(x) = ^V'j(x). Then, 



M 



< cN a/d , / f(x)dx = 1. 
JM 

The estimate of the W (M) norms of (J + A) a/2 ipj(x) and (I + A) a/2 f(x) when a/2 
is an integer follows from the fact that (/ + A) a ^ 2 is a differential operator and the terms 
(/ + A) Q//2 ipj(x) have disjoint supports. When a/2 is not an integer the estimate follows 
by complex interpolation. Anyhow, the case p = 2 is elementary. If 5 > 1 and a5 is an 
integer, 

1/2 

£(l + A 2 )V/(A)f ' 



< 



x (l-l/S)/2 , n 1/25 



□ 
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3 Further results 



The following is Result (5) in the Introduction and it states that a quadrature rule which 
gives an optimal error in the Sobolev space W a ' 2 (/A) is also optimal in all spaces W^' 2 (A4) 
with d/2 < (3 < a. 

Theorem 3.1. If dv[x) is a probability measure on M., then the norm of the measure 
dv(x) — dx as a linear functional on W a,<2 (.M) decreases as a increases. Moreover, if 
the norm of dv[x) — dx on W a ' 2 {M) is r~ a , 



f(x)dv(x) — / f(x)dx 
m Jm 



< r 



then for every d/2 < (3 < a there exists a constant c which may depend on a, (3, M., but 
is independent of r and dv[x), such that 



f[x)dv{x) 



f(x)dx 



Proof. By Corollary l2.4l (2) the norm of the measure dv(x) — dx as a linear functional on 

W a - 2 (M) is 



B la (x,y) dv(x)dv{y) - 1 



M JM 



) 1/2 = {£(i + A 2 n^(A)r 

^ I A>0 



1/2 



Since (1 + A 2 ) a < (1 + A 2 ) p when < a, it follows that this norm is a decreasing 
function of a. Write du(x) — dx = dfi(x). By Lemma [2T6l (D. the norm of the functional 
Im on W°' 2 (Ai) can be written as 



B a (x, y) d/j,(x)dfi(y) 



M JM 



1/2 



\v(aY l f r- 1 exp (-t) ( [ [ W(t,x,y)dii(x)dp J (y))dt 
I Jo \Jm Jm J 



1/2 



Assuming that this norm is r~ a , one has to show that the corresponding expression with 
f3 instead of a is at most cr~^ . Since (3 < a, the integral over 1 < t < +oo satisfies the 
estimate 



r -1 exp(-i) / / W (t,x,y)dfi(x)d^(y) ) dt 
i \Jm Jm J 

< [ t a_1 exp(-t) ( f [ W(t,x,y)dfi(x)dfj,(y)) dt 
Ji \Jm Jm J 

<T(a)r~ 2a . 
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Similarly, since /3 < a the integral over r 2 < t < 1 satisfies the estimate 

t^' 1 exp (-£) ( f f W{t,x, y) dfj,(x)dfj,(y)) dt 

1 \JM JM J 

<r 2a ~ 2fi [ t"" 1 exp (-t) ( / / W(t,x,y)dft{x)dfi(y) J 



Finally, by the Gaussian estimate on the heat kernel in the proof of Lemma l23l if < t < 
r~ 2 then 

t d/2 W (t, x, y) < cr~ d W (r~ 2 , x, y) . 
It then follows that if /3 > d/2 the integral over < t < r~ 2 satisfies the estimate 



t 13 ' 1 exp (-*) / / W{t,x,y)dfj,{x)dfj,{y) ) dt 
<cr~ w I W (r~ 2 ,x,y) d\y\ (x)d\fi\ (y). 

JM JM 

It remains to show that the last double integral is uniformly bounded in r. Since d \fx\ (x) = 
dv(x) + dx and since f M W (r~ 2 , x, y) dx = 1, replacing d \n\ (x) with du.{x) it suffices 
to show that 

W (r~ 2 ,x,y) du,(x)du,(y) < c. 



' M JM 

By the assumption on dfx(x) and the eigenfunction expansion of W (r~ 2 , x, y), 

W (r~ 2 ,x,y) d^{x)dn{y) 

M JM 



< r 



W (r 2 ,x,y) dfx(y) 



M 



■ 1/2 
|2 ' 



< |E (1 + l^(A)| 2 | 7 sup { (1 + A 2 )" exp (- (A/r) 2 /2) } . 

Finally, the last sum with {.F/^A)} is the norm of the measure dfx(x) as functional on 
W a ' 2 (/A), hence by assumption it is r~ a , and the last supremum is dominated by r 2a . □ 

As we said, the above result is only one way, from a to /3 < a. If the norm of 
du(x) — dx on W a ' p {M) is r~ a and if f3 > a, then one cannot conclude that the norm 
of du(x) — dx on W^' p (Ai) is at most cr _/3 . As a counterexample, it suffices to per- 
turb a good quadrature rule with nodes {zj}^ =1 and weights {ouj}j =1 by moving the last 
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point zn into a new point t^, so that the new quadrature differs from the old one by the 
quantity uj n \ f (z N ) — f (t N )\. If a > d/p + 1 then the function / is differentiable and 



*/un one ob- 



un \f (z N ) - f (*jv)I » w ^ Isw - *jv| • Then, by choosing |zjv 
tains a quadrature rule which gives an error m r~ a in all spaces p (/A) with (3 > a. 
The counterexample when d/p < a < d/p + 1 is slightly more complicated but similar. 

In all the above results, the accuracy in a quadrature rule has been estimated in terms of 
the energy of a measure. It is also possible to estimate this accuracy in terms of a geomet- 
ric discrepancy. The Bessel kernel can be decomposed as superposition of characteristic 
functions, 



B a (x,y) 



X{B a (x,y)>t}(x)dt. 



If 1 < p,q < +oo and 1/p + 1/q = 1, by Corollary 12 .41 and Minkowski inequality, the 
following Koksma Hlawka type inequality holds: 

1/9 



f(x)djj(x) 



M 



< 



M 



B a (x, y)dfi(x) 



M 



< 



M 



M 



X{B"(x,y)>t}{x)d/l(x) 



dy 



dy 



dt. 



The quantity \ j M X{B a (x, y )>t}(x)dfi(x)\ is the discrepancy of the measure dfi(x) with 
respect to the level sets {B a (x,y) > t}. It can be proved that, for specific measures and 
at least in the range < a < 1, the above estimates are sharp and they can lead to optimal 
quadrature rules. In particular, the following is Result (6) in the Introduction. 

Theorem 3.2. Denote by 5{t) the diameter of the level sets of the Bessel kernel {B a (x, y) > 
and assume that there exists r > 1 such that the discrepancy of the measure dp,(x) with 
respect to {B a (x, y) > t} satisfies the estimates 



M 



X{B<*(x, y )>t}{x)dn{x) 



< 



r- d if5(t)<l/r, 
r- l 5{t) d - 1 if8(t) > l/r. 



Also assume that 1 < p < +oo and a > d/p. Then there exists a constant c, which may 
depend on a and p and on the total variation of the measure \p\ (M), but is independent 
ofr, such that 



f(x)dp(x) 



M 



cr ~ Q ||./llw^ 5f0<a<l, 
<< cr- 1 log(l + r) ||/|| w „,, if a 

cr ~ l ll/llw°.p l f a > 1 - 



I, 



Proof. Observe that the above hypotheses on the discrepancy match the estimates in 
Lemma [2791 Indeed, by this lemma, the measures du(x) which give exact quadrature for 
eigenfunctions with eigenvalues A 2 < r 2 have discrepancy 

~ d ifs<l/r, 



ii 



dvix) — 



{\x-y\<s} 



dx 



{\x-y\<s} 



< 



cr 



cr 1 s d 1 



if s > l/r. 
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Actually, these estimates hold not only for balls {\x — y\ < s}, but also for sets with 
boundaries with finite d — 1 dimensional Minkowski measure, such as the level sets 
{B a (x,y) > t}. Also observe that these estimates are natural, since the discrepancy of 
large sets is qualitatively different from the one of small sets. If 1 < p, q < +oo and 
1/p + l/q = 1, by Corollary 12 .41 and Minkowski inequality, 



f(x)dfi(x) 



M 



< 



M 



M 



X{B<*{x,y)>i}{x)dn{x) 



1/9 



dy \ dt. 



By Lemma |231 when < a < d then B a (x, y) \x — y\ a ~ d , the level sets {B a (x, y) > t} 
have diameters 5(t) w min {l, tVt"-* 4 ) | and the boundaries {B a (x,y) = t} have surface 
measure of the order of <5(t) d_1 ~ min |l, t( d_1 )/( a_d )|. Hence the estimate of the dis- 
crepancy of small level sets with t > r d ~ a gives 



l M 



M 



X{B"{x,y)>t}(x)dfi(x] 



dy 



< sup 

y£M 



' M 

< sup 

y£M 



X{B"ix,y)>t}{x)djJ,(x) 



(9-l)/9 



M 



X{B<*(x,y)>t}(x)d[l(x) 



Jm Jm 

(9-i)/a 



X{B«(x,y)>t}(x)d\ii\ (x)dy 



1/9 



{c\fi\(M)t d ^} 1/q 



< cr -%- 1 )/9^/9(a-d) 



Hence, if a > <i/p the integral over r a < t < +oo satisfies the inequality 



M 



X{B"(x,y)>t}{x)dll(x) 



M 



dy > dt 



<cr -d(q-i)/q / t d/q{a - d) dt < cr~ a . 



Similarly, the integral over < t < r d a , that is the discrepancy of large level sets, 
satisfies the inequality 



X{B"(x,y)>t}(x)dfi(x) 



1/9 

dy dt 



< r 



o UM jm 
•d-a ( cr- a if < a < 1, 

min {l,t {d - 1)/{a - d) } dt < < cr- 1 log(l + r) ifa = l, 

I if a > 1. 



C:?' 



The proof in the case a = d is similar and it follows from the estimate B a (x,y) ~ 
— log (\x — y\). The proof in the case a > d is even simpler, since in this case B a (x, y) 
is bounded and it suffices to integrate on < t < sup x eM B a (x, y) the inequality 

\S M X{B"(x, y )>t}(x)dfj,(x)\ < cr-\ □ 
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In particular, it follows from Lemma [279l Theorem 12. 12[ Theorem 12. 141 that, at least 
in the range < a < 1, Theorem |3. II gives an optimal quadrature. We conclude with a 
series of remarks. 



Remark 3.3. As we said, the assumption a > d/2 withp = 2 in Theorem 12 .7 [ or a > d/p 
with 1 < p < +oo in Theorem 12.121 guarantees the boundedness and continuity of 
/ (x), otherwise / (zj) may be not defined. This follows from the Sobolev imbedding 
theorem. Indeed, the imbedding is an easy corollary of Lemma l2.6l A function f(x) is in 
the Sobolev space W a,p (/A) if and only if there exists a function g{x) in LP (M) with 

f(x) = [ B a (x,y)g(y)dy. 

JM 



When 1 < p, q < +oo, l/p + 1/q — 1, d/p < a < d, then B a (x, y) < c\x — y\ a ~ d is in 
L 9 {M) and this implies that distributions in the Sobolev space W a ' p {M) with a > d/p 
are continuous functions. Indeed they are also Holder continuous of order a — d/p. 

Remark 3.4. When the manifold is a Lie group or a homogeneous space, one can restate 
Theorem [2J] in terms of convolutions. In the particular case of the torus T d = R d /Z d , let 

A{x) = ip{k) exp (2,7rikx) , B{x) = ip{k)~ l exp (2iiikx) . 

k&z d 



Then, if 1 < p,q,r < +oo with 1/p + 1/q = 1/r + 1, 

r > 1/r /• J, 

l/i(y) dx\ = \ 
<<l \A*f{x)\ p dx\ ' { [ \B*fi(x)\ q dx 



B * A* f * /i(x)| r dx 



1/9 



1/r 



In the case of the sphere E> d = {x E M. d+1 , \x\ = l}, let {Z n (xy)} be the system of 
zonal spherical harmonics polynomials and let 



M X V) = ^2^P( n ) Z n (xy) , B(xy) = ^2i>(n) 1 Z n (xy) . 



n=0 



n=0 



Then, if 1 < p,q < +oo with 1/p + 1/q = 1, 



/ 



f (x) dfi(x) 



< 



A( x y)f(y) d y 



dx 



B(xy)d/i{y) 



dx 



1/Q 



Both results on the torus and the sphere follow from Young inequality for convolutions. 
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Remark 3.5. A result related to Theorem 12.11 is the following. Identify T d with the 
unit cube {0 < Xj < 1} and denote by XP(y)( x ) me characteristic function of the par- 
allelepiped P{y) = {0 < Xj < Uj}. Then define 



P d 

B(x)= / X P(y)(x)dy-2- d = l[(l 



Xj) - 2 



fcez d -{o} \ \kj=o J \fc^o 
Also, define the differential integral operator 



exp (2irikx) . 



a * f(x) = n 2 n 2mk i ex p ( 27r ^ x ) 

k^O \kj=0 J \kj^0 J 

KiXd^ 1 J i±j l<i^i<rf J Mi, 7 



9 d 



■fix)- 



dx\...dxd' 

Observe that, as in Theorem 12. 1[ the Fourier coefficients of the distribution A(x) and 
of the function B(x) are one inverse to the other, however here the Fourier coefficients 
are function of the lattice points 2mk, and not of the eigenvalues An 2 |/c| 2 . If dv(x) = 
N^ 1 J2f=i ( x )> an d ^ 1 — Pi Qi r — +°° w i m 1/p + 1/? = 1/t + 1, then 

r s 1/V 



{/ N-^f{x-Zj)-j 
yjf d j =l Jj d 



f{y)dy 



dx 



< I I \A*f(x)\ p dx 



i/p 



\B * u(x)\ 9 dx 



1/9 



The norm of A * f(x) is dominated by an analogue of the Hardy Krause variation, 



\A*f(x)\ p dx 



i/p 



i<j<d l/ T 

+2 d - 2 Yl If 







i dx j 
jd-2 dxidxj 



f( x ) \\dx r 



i= 3 



dxj 



-f(x) Y[ dx h 



fd 



dxi...dx c 



dx 



1/p 



26 



L. Brandolini, C. Choirat, L. Colzani, G. Gigante, R. Seri, G. Travaglini 



The norm of B * u(x) is dominated by the discrepancy of the points {zj} =1 with respect 
to the family of boxes P(y), 

\B * u(x)\ q dx > 
jd J 

N d 

N ~* Xp (y) ( z i + x ) - II yj 

3=1 3=1 

In particular, the case p — 1 and g = +oo is an analogue of the Koksma Hlawka inequal- 
ity. See [|20ll . A generalization of this classical inequality is contained in J6J. 





Remark 3.6. By Lemma [2T6l (l). the Bessel kernel B a (x, y) with a > is a superposition 
of heat kernels W (t, x, y). Indeed, it is possible to state an analogue of Corollary 12.41 in 



terms of the heat kernel, without explicit mention of Bessel potentials: If {zj} N =l is a 



sequence of points in Ai, if - =1 are positive weights with J2j u j = 1» an d if f{x) is 
a function in W Q,P (A^) with a > d/2, then 



< <^ r (a) 



+oo 



- r 

j=1 Jm 



N N 
4=1 j=l 



f(x)dx 



1/2 



t a-1 exp (-t) dt 



This suggests the following heuristic interpretation: Mathematically, a set of points on 
a manifold is well-distributed if the associated Riemann sums are close to the integrals. 
Physically, a set of points is well-distributed if the heat, initially concentrated on them, in 
a short time diffuses uniformly across the manifold. 



Remark 3.7. In order to minimize the errors in the numerical integration in Corollary l2.4l 
(3), one has to minimize the energies 



N N 



/ / B 2a (x,y)dv(x)dis(y), VVw^S 
Jm Jm i=1 j=l 

These are analogous to the energy integrals in potential theory 



[Zi, Zj) . 



\x — y\ du(x)du(y). 

I M JM 

See [14J. When d < a < d + 1 the kernel B 2a (x,y) is positive and bounded, with 
a maximum at x = y and a spike A — B \x — y\ 2a ~ d when x — > y. In particular, the 
gradient at x = y is infinite. This implies that in order to minimize the discrete en- 
ergy J2i j UiUjB 2a (zi, Zj) the points {zj} have to be well separated. This suggests the 
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following heuristic interpretation: Mathematically, a set of points on a manifold is well- 
distributed if the energy is minimal. Physically, a set of points, free to move and repelling 
each other according to some law, is well-distributed when they reach an equilibrium. 

Remark 3.8. It can be proved that if 2a > d + 2 then 

\B 2a (x,x) -B 2a (x,y)\ < c\x-y\ 2 . 

This estimate in the proof of Theorem [277] yields that for most choices of sampling points 

Zj E Uj, 



N . 

ujjf (zj) — J f{x)dx < c ^ a ^- N jdiameter (U 



sd/2+1 



The same result holds if 2a = d + 2, with a logarithmic transgression. Observe that 
these estimates hold for most choices of sampling points, but not for all choices. Indeed, 
if the manifold Ai is decomposed in disjoint pieces Ai = U\ U U2 U ... U Un with 
measure aN- 1 < \Uj\ = uj < bN- 1 and diameter (Uj) < cN' 1 ^, if f(x) is a smooth 
non constant function, and if the points Zj E Uj are the maxima of f(x) in Uj, then 



Y^jLi Ujf ( z j) ls an upper sum of the integral J M f(x)dx and 

- r 

J2cjjf( Zj )- / f(x)dx > ciV- 1 ^. 

Remark 3.9. Theorem 13.21 gives an estimate of the accuracy in a quadrature rule in 
terms of the discrepancy of a measure with respect to level sets of the Bessel kernel. The 
following argument shows that when the manifold is a sphere, or a rank one compact 
symmetric space, then the level sets of the heat kernel {W (t, x, y) > s}, and hence of the 
Bessel kernels {B a (x, y) < t}, are geodesic balls {\x — y\ < r}. The Laplace operator 
on the sphere E> d with respect to a system of polar coordinates x = a), with < $ < n 
the colatitude with respect to a given pole and a E S d_1 the longitude, is 

^ = A (lV) = - sin 1 ^ (tf) A ( sin ^ ^ + K . 

Let u (t, x) be the solution of the Cauchy problem for the heat equation 

d 

—u(t,x) = -A x u (t,x) , 
u(0,x) = f(x) . 

If / (x) depends only on the colatitude i3, if it is an even function decreasing in < i3 < n, 
then also u (t, x) depends only on the colatitude and it is an even function decreas- 
ing in < •§ < 7r. In order to prove this, set u(t, x) = U(t, f (x) = F ($), and 
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sm d ' 1 (0) dU(t, 0)/<90 = V(t, 0). Then 



d0 



[7(0, 0) 



d_ 



0_ 

Of 



92 



,d-l 



d 



sin(0) d-d 



V{0, 0) = sin' 1 - 1 (0) ^F(0), 
[ 7(t,0) = y(t,7r) 



90' 
0. 



If F (0) is decreasing in < < it, then V(0, 0) < and, by the maximum principle, 
V(t,'&) < 0, hence U(t, 0) is decreasing in < < n. In particular, by considering 
a sequence of initial data {f n (x)} which depend only on the colatitude 0, even and de- 
creasing in < < it, and which converge to the Dirac 5(x), one proves that the heat 
kernel W (t, cos (0)) is decreasing in < < it. Since Bessel kernels are superposition 
of heat kernels, they are also superposition of spherical caps. 

Remark 3.10. In Q and [|22ll the discrepancy of orbits of discrete subgroups of rotations 
of a sphere are studied. Let Q be a compact Lie group, /C a closed subgroup, M. = Q/K, a 
homogeneous space of dimension d. Also, let H be a finitely generated free subgroup in 
Q and assume that the action of % on M. is free. Given a positive integer n, let {(?j}^ =l 
be an ordering of the elements in % with length at most n and for every function f{x) on 
M, define 

N 

Tf(x) = N- 1 J2f(^)- 

This operator is self-adjoint and it has eigenvalues and eigenfunctions in h 2 (A4). More- 
over, since the operators T and A commute, they have a common orthonormal system 
of eigenfunctions, A(p x (x) = \ 2 (f\(x) and Ty? A (x) = T(\)ip\(x). All eigenvalues of T 
have modulus at most 1 and indeed 1 is an eigenvalue and the constants are eigenfunc- 
tions. Assume that all non constant eigenfunctions have eigenvalues much smaller than 1. 
Then, if a > d/2, 



N 



OjX) 



J'=l 



/ f(x)dx 




1 M 





A^O 
1/2 



< {sup {|T(A)i>} jx: (i + xr i^/wi 2 } 1 |e (i + at q i^(x)i 2 | 7 



<c<jsu P {|T(A)|} 

A^O 



(I + A) a/2 f(x)\ dx 



1/2 
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The absolute convergence of the above series is consequence of the Sobolev's imbeddings, 
or the Weyl's estimates for eigenfunctions. In particular, when M. = SO(3) / SO(2) is 
the two dimensional sphere and % is the free group generated by rotations of angles 
arccos(— 3/5) around orthogonal axes, it has been proved in G2l that the eigenvalues of 
the operator T satisfy the Ramanujan bounds 



sup{|T(A)|} < cN- 1 ' 2 \og{N). 

A^O 



Hence, for the sphere, 



N 

, =1 Jm 



f(x)dx 



< cN~ 1/2 log(iV) 



M 



(I + A)"/ 2 f{x 



dx 



1/2 



All of this is essentially contained in Il22l . Although this bound iV -1 / 2 log(iV) is worse 
than the bound N~ a ^ 2 in Corollary 12.131 the matrices {aj} have rational entries and the 
sampling points {<?jx} are completely explicit. 
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